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We analyze analytically and numerically the origin of the singularity in the course of the collapse 
of a wormhole with the exotic scalar field with negative energy density, and with this field 
together with the ordered magnetic field H. We do this under the simplifying assumptions of the 
spherical symmetry and that in the vicinity of the singularity the solution of the Einstein equations 
depends only on one coordinate (the homogeneous approximation). In the framework of these 
assumptions we found the principal difference between the case of the collapse of the ordinary scalar 
field $ with the positive energy density together with an ordered magnetic field H and the collapse 
of the exotic scalar field <& together with the magnetic field H. The later case is important for the 
possible astrophysical manifestation of the wormholes. 
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I. INTRODUCTION 

Wormholes (WHs) are hypothetical short topological 
tunnels connecting two different distant asymptotically 
flat regions of the Universe, or such regions belonging 
to different universes in the model of the Multiverse [l[. 
They are typical relativistic objects. 

The problem of wormholes has a long history. In the 
framework of General Relativity, first attempts to con- 
struct such objects were performed by Flamm and 
Einstein and Rozen [3]. Later on, many theoretical as- 
pects of the WH problem were investigated, see , for ex 
ample EL [EJEH ijjlpjlljll 

19, 2a 21, pi pi pi HE HE S3, H 

3a . l3q . [371 [3a . KM 140 ] . During the last few decades, inter 
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est in wormholes increased in connection with the early 
Universe [ill l42l . |4 3|. [43 an d various problems in physics 
and astrophysics [J EE EE S3, El IS HI EL El EE H • 

In the works [IE EL Q3 the hypothesis that some 
known astrophysical objects (e.g. quasars and active nu- 
clei of some galaxies) could be entrances to wormholes 
was considered. Wormholes may have existed as pri- 
mordial objects in early stage of the expanding Universe 
8, 43]. It is possible that such primordial wormholes 
could be preserved after the end of the inflation [H EI • 
This hypothesis can explain some observable facts in as- 
trophysics and can predict new phenomena [ll| ■ 

In the middle of the last century it was shown in Gen- 
eral Relativity that a vacuum WH pinches so quickly 
that it cannot be traversed even by a test signal mov- 
ing with the velocity of light (see review in fl5|). In 
order to prevent the shrinking of a WH and to make 
it traversable, it is necessary to thread its throat with 
so-called exotic matter which is matter that violates the 
averaged null energy conditions (see [E [lfj, [lg, [lTj ) . Dif- 
ferent types of wormholes may exist depending on the 
type of exotic matter in their throats [H HE El . For 



example it could be a "magnetic exotic matter" in which 
the main component is a strong ordered magnetic field 
plus a some amount of a "true exotic matter" [ll[. An- 
other type is a "scalar exotic matter" in the form of a 
scalar field with a negative energy density [E HE] ■ One 
more type is a mixture; "magnetic-negative dust exotic 
matter" which is a mixture of an ordered magnetic field 
and dust (matter with zero pressure) with negative mat- 
ter density [2l| . The physical properties of different types 
of WHs are different. WHs with scalar exotic matter 
were the subject of very intense investigations both an- 
alytically and numerically. Possible dynamics of such 
WHs has been analyzed analytically for example in pa- 
pers P, HE Hi HE L2E H3, HE HI, lafll ■ The most impor- 
tant, however, was the numerical analysis [H, HE HE EH- 
It was shown that the static WHs of this type are unsta- 
ble. Perturbations trigger the evolution of the WH and 
the evolution of the exotic scalar field which maintains 
it. As a result, the WH either collapses or expands. In 
the case of the collapse - a black hole (BH) arises. If 
the collapsing WH has both the exotic scalar field and 
the magnetic field, the structure of the singularity inside 
the resulting BH requires special investigation. The rea- 
son for extra attention is that the magnetic field in WHs 
or their remnants has special manifestation in the astro- 
physical observations in the hypothesis of the existence 
of WHs in the Universe [ll| • 

The goal of this paper is to analyze (under simplifying 
assumptions) the nature of the singularity created during 
the collapse of the WH with the exotic scalar field and 
with the magnetic field. 

The paper is organized as follows. In section |TT] we 
describe a model in the framework of which we analyze 
the problem. In the subsection III Al the equations are 
written out. In the subsection III Bl we use the leading 
order analysis for the case of a WH without magnetic 
field. In the subsection III CI we use a numerical code to 
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solve the equations to understand the behavior of the 
model for several sets of its parameters. In subsection 
III Dl we use the leading order analysis for the case with 
the magnetic field. In the subsection III El the numerical 
analysis is applied to the case with the magnetic field. In 
the section HTT1 we discuss the results. 
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II. 



THE MODEL 



We consider collapse of a WH with the formation of 
a BH. The main feature inside a BH is its space-time 
singularity. Our goal is to investigate the nature of this 
singularity, using some simplifying assumption. First, 
we consider the spherical WH and BH nonlinearly per- 
turbed by a minimally coupled and self-gravitating mass- 
less exotic scalar field $ with the negative energy den- 
sity e < 0. We will also consider the same case but with 
an additional radial magnetic field. It was shown (see 
[EE IE3, HH) that, in the considered model, in the close 
vicinity of a space-like singularity of a BH, all processes, 
as a rule, have high temporal gradient (much higher than 
the spatial gradients along the singularity) and that the 
processes depend on the properties of a very restricted 
space region. It follows from here that, for clarification 
of some physical processes, one can use a homogenous 
approximation and that all processes and geometry de- 
pend on the time coordinate only. We assume also that 
in the close vicinity of a time-like singularity of a BH 
one can use a homogenous approximation also, but now 
all processes and geometry depend on the radial space 
coordinate only. 



where the primes denote differentiation with respect to 
r. Tensor T represent here contribution from the $ field 
and tensor H represents contribution from a free radial 
magnetic field. In the WH, H is a sourceless magnetic 
field captured by the topological structure of the 3-D 
geometry. For both the WH and arising from the collapse 
of it BH the components of the tensor H are given by: 
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where the constant q characterizes the strength of the 
magnetic field. The exotic scalar field \& is governed by 
the Klein-Gordon equation \I>:^ = (; denotes the covari- 
ant derivative), whose first integral reads 



-g rr sin © 



(7) 



where d is a constant and g is the metric determinant. 
For the exotic field, the value of d is pure imaginary, 
d 2 < 0. Of course <&' does not depend on O because the 
metric determinant g has a factor sin 2 O . 

The components of the tensor T for the exotic scalar 
filed * are (see [Hj]): 



Tl = -e 



(8) 



A. The equations 

We start with the general homogeneous spherically 
symmetric line element: 



ds l 



g u {r)dt 2 + g rr (r)dr 2 + r 2 dQ? 
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Inside a BH the region between the event horizon (EH) 
and the Cauchi horizon (which exists in the case with 
the radial magnetic field; we will call it horizon-2 (H-2)) 
is so-called T-region (see [l5[ and references therein). In 
this region r is time- like and t is space- like coordinate. To 
describe the contraction, we should consider the variation 
of the time coordinate r from bigger to smaller values. 
The r — r, t — t and — components of the Einstein 
equations (with c = 1, G = 1) are given by (see (58[) 
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where eo, gtt,init, finite d 2 < are all constants. Sub- 
stitution of fo] 151- [T2f into (0]), (gD enables us to find 
the unknown functions g rr {r) and <?tt(r). Equation ([5]) is 
the consequence of ([3]), (UJ) and hence can be used as a 
control of the calculations. 



gtt - grrgtt + rgi u 
r 2 g rr gtt 



8tt(T; + H r r ) 



(3) 



B. Leading order analysis for the case q — 



We will use the method which has been proposed by 
Burko [Hf! [H?! and after that used in [58|. This is the 
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non-linear generalization in the homogenous case of lin- 
ear analysis of In this section we consider the case 
<7 = 0. Let us consider the leading order terms in a se- 
ries expansion for the metric functions and the leading 
order terms in the Einstein equations ([3]), (J3J, q — near 
the singularity r = 0. We assume that the leading terms 
for each g tt and g rr can be written in the power form 
(const ■ r P ower y We demand also that the correspond- 
ing terms in the Einstein equations (O, ((4]) tend to zero 
when r — ► 0, otherwise this method doesn't work (see 
Appendix [X}. In the vicinity of r = the solution of 
([5]), ([4"|1 can be written, as the first approximation, in the 
following form 



^ = -0 + 2)^, 



(13) 



(14) 



where m, C and /? are all constants. For the scalar field 
^, the leading order (in r) analysis gives: 



(i) 

M 



V/5+Tlnr. 



We have also for the constant (i 2 (see ([7]) and (|12p) 



09 + 2) 



(15) 
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Here we use constants mo and C which were introduced 
in (58| . These constants have direct physical meaning in 
the case of empty (without \& field) BH, when /3 = — 1. 
In such a case, m is the BH mass, C is a gauge param- 
eter related to the possibility of changing the scale of 
measurement of the t space coordinate. 

In the case of presence of the ^ field, (3 determines the 
strength of this field. In contrast to the case analyzed by 
Burko in fHH, , in our case here ^ should be imaginary 
and from (11611 we have the restriction: 



(3 < -1. 



(17) 



On the other hand, for too > we obtain from ((T 
and JTH): 



C > 0, 



(3 + 2 > 0. 



(18) 



(19) 



Thus we have for (3 the restrictions — 2 < (3 < — 1. In the 
case of violation of the inequality (JX9J) this method does 
not work as we mentioned above and in Appendix A. 

The above solution (TTB")) , (fT4"|) is generic, because it 
depends on three arbitrary parameters. In our case m 
and (3 are essential physical parameters, and C is a gauge 
parameter. Using the arguments which are analogous to 
the arguments of Burko [561, (57j one may conclude that 
the dependence on two physical and one gauge parameter 



means that the solution is generic. If we suppose that in 
the solution (TTB"]) . (fT4"|) mo < 0, C > 0, then it means 

that <?jj < 0, grr > 0; t is the time coordinate and r 
is the space coordinate. In this case r = is a time-like 
singularity. Probably it corresponds to negative mass of 
the object and a naked singularity. 



C. Numerical analysis 

We use a simple numerical code to solve numerically 
equations ([3]), (UJ) to understand the behavior of model 
(H|) , ^ as a function of the parameters of the model. 

We will represent the results for the case of the exotic 
scalar field ^ with e < 0, and compare them with the 
results for the case of the ordinary scalar field with 
e > 0. For the latter case <fr-field with e > 0, we consider 
also the case of the presence of the magnetic field g ^ 
to emphasize the essential difference the cases q = and 

We start the computation from ri n u — 0.95r\E#, and 
set qtt,init) qrr,init equal to their values at for the zero 
matter content Schwarzschild (or Reisner-Nordstrom if 
q =/= 0) solution with initial mo = 1 and q — 0.95 if 
q y£ 0. Here teh is the value of the event horizon for the 
zero matter content BH with the same parameters. We 
consider several values of eo, which is the characteristic 
of the initial amplitude of the scalar field. 

The parameters of the model, including the exponent 
(3, depend on the value of the scalar field eo- We start 
from the extremely small amplitude |eo| and consider the 
cases (see Figure [1]): 



£ = +0.0001, q = 0, 
eo = -0.0001, q = 0, 
e = +0.0001, q = 0.95, 



(20) 
(21) 
(22) 



In Figure 1(a) one can see the propagation (r vs. t) of 
the incoming signal with the velocity c (c = 1). Figure 



1(b) shows the mass function (see |58j): 



m 



— 



fJr 



(23) 



Figure 



1(c) presents the evolution of the metric functions, 
and Figure 1 (d) shows the evolution of the 
(3 and d\og\g rr \/d\ogr_=_a. 



g tt and g rr , 
values d log gtt /d log r 

It is clear from the results shown in Figures 1(c) - 1(d) 
that for this extremely small amplitude |eo] m the 
case q = the propagation of the light-like sig- 
nal, evolution of gtt and g rr and dlogg t t/dlogr = (3 
and dlog\g rr \/dlogr = a are practically the same for 
eo = ±0.0001 and the difference is not visible in the plots. 

On the other hand, the behavior of the mass function 



for the cases e = 0.0001 and e = -0.0001 (Figure [1 (a) [ ) 
is quite different. Of course the reason for this is the 
opposite signs of e for these two cases. 

The physical processes which can lead to a nonlinear 
change of the mass function are the following: 
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(c) Metric functions gtt and |g r r| vs r. 

FIG. 1: The cases of e = 0.0001, q = 0, 
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(b) Mass function versus r. 
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(d) The d log gtt /d log r and dlog\g rr \/dlogr vs 
-0.0001, q = 0, and e = 0.0001, q = 0.95. 
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A) The mass m inside a sphere can change because of 
the work of pressure forces on the surface of the 
sphere, 

B) The mass can change due to the mass inflation (58|. 

Both processes are described in [58[- The process B) 
is important near the horizon-2 (H-2) in the case of pres- 
ence of the magnetic field q ^ and e > 0. In the case 
q = horizon-2 (H-2) does not exist. Thus in the case 
q = the process B) is not important, and variations 
of mass function is small and opposite for the two cases 
e = 0.0001 and e = -0.0001. 

The case of eq = 0.0001, q = 0.95 is quite different. 



Figure 1(a) shows that in this case the light- like signal 
goes close to the initial value of r for the (H-2) horizon 
r = rn-2 during a long period, and only after that it 
comes to r ~ (beyond the framework of the figure). 
Note that the light-like signal does not cross the H-2 
horizon. The horizon itself shrinks down to r = un- 
der the focusing effect, as it is described in [58]. The 
metric functions gtt and g rr go to zero very fast, and 
very large asymptotic values d log g tt /d log r = (3 and 



d\og\g rr \/d\ogr — a are clearly seen in Figure 1(d) 
When the solution is close to r — th-2 the mass infla- 



tion manifests itself (see Figure 1(b) ), the mass function 
increases dramatically. 

We next study the solutions for the larger values of | eo I ■ 
We start with the case eo > and consider the following 
values (see Figure [2]) : 



e = 0.011, q = 0, 
e = 0.012, q = 0, 
e = 0.050, q = 0, 



(24) 
(25) 
(26) 



The propagation of the light-like signal (Figure 2(a) ) 
qualitatively is the same for all cases but for larger val- 
ues of eo it comes to r = later, behavior of the metric 



functions, Figures 2(b) 2(c) is interesting. The analysis 



of the leading order terms of the series expansion for the 
case e > (this means that d 2 > 0) leads to the conclu- 
sion that (3 > — 1, but [3 can be positive and negative. 
Numerical analysis shows that for rather small eo < 0.011 
the function g u increases when r — > 0, but for the larger 
eo > 0.012 the function |g tt | decreases when r — > (see 
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(a) r versus t. 



(b) Metric functions gtt and \g rr \ vs r. 
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(c) The d log gtt/d log r and dlog \g rr \/dlog r vs r. 



(d) Mass function versus r. 



FIG. 2: The cases of e = 0.011, q = 0, e = 0.012, q = 0, and e = 0.050, q = 0. 



TABLE I: The limit 
(dlog \g rr \/d log r) = c 



values of (dlog g t t /dlog r) 
: for q — 0, for r — > 0. 



/? and 



eo 


a 


P 


0.0001 


1.00 


-1.00 


0.0010 


1.01 


-0.99 


0.0025 


1.02 


-0.98 


0.0005 


1.06 


-0.94 


0.0100 


1.52 


-0.48 


0.0110 


1.79 


-0.21 


0.0120 


2.16 


+0.16 


0.0200 


8.95 


+6.95 


0.0500 


53.8 


+51.8 



also Tabic HJ. 

As it can be seen from Table U the function lim r ^o P 
where (3 = d log gtt /d log r increases monotonically when 
6q varies from f3 = — 1 when eo — 0, to f3 = +51.8 when 
eo = 0.05. The critical value of eo, when (3 changes 
its sign, is eo. cr it ~ 0.0115. The value lim r ^ ot where 
a = d\og \g rr \/dlogr also increases monotonically when 
eo varies from a — 1 when eo = to a — 53.8 when 



eo = 0.05. The sign of a never changes. 

The difference in behavior of the metric functions leads 
to difference in the properties of the mass function. As we 
mentioned above, in the case q = the horizon H-2 does 
not exist and the mass inflation (the process B)) is not 
important for the evolution of the mass function. Thus, 
the process A) , which is related to the deformation of the 
volume of the reference frame, plays the main role. The 
longitude deformation is proportional to (gtt) 1 ^ 2 , while 
the transversal deformations are proportional to r 2 . In 
all cases it leads to the increase of the mass with decrease 
of r. For large positive [3 corresponding to large eo this 
increase is very fast. 

We now turn to the discussion of the case e < 0, q = 
(see Figure [3]), which corresponds to the case of the col- 
lapse of the WHs analyzed numerically in (55|. Here 
in the limit r — > 0, j3 is always less than -1. With 
the increase eo the limit lim r ^o \P\ increases insignifi- 
cantly. Dropping the absolute value, the limit varies from 
(3 = -1.00 at e = -0.0001 to j3 = -1.01 at e = -0.05. 
The behavior of the metric functions g u , g r r is similar 
for all cases 0.0001 > e > —0.05, see Figures 3(a 



3(b) 
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The light-like signal comes to r = later for the larger 
values of |eo|, see Figure 3(c) This property is opposite 



to the case of eo > 0. The changes of the mass function 
m for the cases eo = —0.01 and eo = —0.05 are also 



insignificant (see Figure 3(d) I, in sharp contrast with the 
case eo > 0. 



D. Leading order analysis for the case q 0. 

We consider the leading order terms in the series ex- 
pansion for the metric functions and the leading order 
terms in the Einstein equations ([3]), (d]), q ^ near the 
singularity r = 0. An assumption of (3 > makes the 
term with the magnetic field q 2 to be negligible, and the 
analysis near r — > becomes the same as for the cases 
considered in the section Hi Bl However, one of our con- 
clusions in section IIIBI was that (3 < — 1 (TIT]) . Thus 
our current assumption of (3 > contradicts the result 
expressed in (fT7|) . The case (3 = also leads to contra- 
dictions. We come to a conclusion that the only allowed 
values for (3 are negative values, (3 < 0, which means that 



in the vicinity of r = we can neglect the scalar field and 
consider d 2 =0, q 2 ^ 0. 

To the leading order (J3J, (U|), d = assume the forms: 



gt rr r - g rr - (g 2 r q 2 )/r 2 = 0, 



(27) 



(28) 



We are looking for solutions to (|27|) . (|28|) in the vicinity 
r = in the form 



g'ttr + gtt + (9rr9ttq 2 )/r 2 = 



} ] = Br? 



flu 



,W = -Ar a . 



Substituting flUD, {H into (|27 |l , ([25 ]) we have 

a = 2, /3 = -2, A= -<T 2 , 



(29) 



(30) 



(31) 



B is arbitrary, but for the correct signature B < 0. 

Of course this is well known the Reissner-Nordstrom 
solution in the vicinity of r = 0. Arbitrariness of the 
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of 




t r 

(c) r versus t. (d) Mass function versus r. 

FIG. 4: The case of e = -0.0001, q = 0.95. 

coefficient B is related to the arbitrariness of the choice FigureS]) and i?-region and come to an internal i?-region? 

of the t coordinate. Note that now g t t < and g rr > 0, To arrive to an answer, one needs to perform a deeper 

this means that t is the time coordinate, r is the space analysis and numerical calculations for the general case 

coordinate, we are in the R- region, see [15|, and r = is when g rr and gtt depend on both r and t coordinates. For 

a time-like singularity. now, we see that here the origin of the r = singularity 

in the T-region is not possible. 

Thus we need to analyze here numerically the case 

E. Numerical analysis e < 0, q 0. when we are in the i?-region and finally come to the time- 
like singularity r = as it was discussed in Section III Dl 

We start with the case q = 0.95, e = -0.0001 and using the leading order analysis. To accomplish that, we 

r init = 0.95r EH . This case differs from those discussed consider the cases (Figure [5J: 
in Section III CI by the sign of e only. The change of sign 
on eo greatly influences the physics of the processes tak- 
ing place. In the case with e > the nonlinear effects e o = -0.0001, q = 0.95, r imt = 0.95r ff _ 2 , (32) 
of mass inflation and focusing of the H-2 horizon down e = —0.0500, q = 0.95, r init = 0.95r#- 2 . (33) 
to r = arized. As a result, both functions \g rr \ and gtt 

finally come to zero at r = 0, and the space- like singular- As we emphasized in Section [II Dl the asymptotic val- 

ity arized (see Figure [T]) in the T-region (see (HI). In our ues of gtt and g rr do not depend on the scalar field. Thus 

current case of e < there are no mass inflation and other the dependence of the solutions on eo is rather weak, 

such effects, and one can observe quite different behavior Note that a and /3 come to their asymptotic values a = 2 

(see Figure^]). The question is, can the collapsing con- and (3 = —2 which do not depend on eo. For larger val- 

figuration cross the border between the T-region (as in ues of | eo I the light-like signal comes to r = at slightly 
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larger t, and the dependence of the mass function m on 
r is stronger (but still rather weak). 

III. CONCLUSIONS 

The evolution of the WHs can lead to their collapse and 
the origin of the singularity. In this paper, we investigate 
the structure of the singularity arising as a result of the 
collapse of the WH with the exotic scalar field and the 
ordered magnetic field. We consider the spherical WHs. 

In the very vicinity of the singularity r — one can use 
an approximation where the solution depends only on one 
coordinate r (the so-called uniform approximation). In 
this region, it is possible to consider the leading order 
terms in a series expansion for the metric functions and 
the leading order terms in the Einstein equations. We 
have demonstrated that in the case of the presence of 
the exotic scalar field only (without a magnetic field) the 
metric functions are 

ff £>~r* g rr ~r^\ g 22 ~r\ (34) 



at r — > 0. Here (3 is a constant 

-2</3<-l. (35) 

We performed the numerical analysis of the uniform ap- 
proximation to understand the behavior of the model as 
a function of its parameters. The results are described in 
SectiorQTBj For the comparison we also investigated the 
case of the collapse with the ordinary scalar field <I> with 
e > instead of the exotic scalar field \& with e < 0. In 
this case instead of inequality f|35[) we have 

0>-l. (36) 

and the numerical analysis shows that f3 changes its sign 
at some critical initial eo i see Table HI 

There is the special interest in the case of the exotic 
scalar field \& together with the magnetic field H. As it 
was mentioned in the Introduction (|T|, the magnetic field 
in the WHs or their remnants has a special manifestation 
in astrophysical observations if the WHs really exist in 
the Universe. 

From the leading order analysis for this case we con- 
cluded that in the vicinity of r = 0, f3 is negative and the 
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magnetic field dominates the scalar field, and the metric 
functions correspond to the Reissner-Nordstrom solution. 
Thus the singularity at r — is in the i?-region and is a 
time-like singularity in the case of our assumptions. Nu- 
merical estimates demonstrate that the asymptotic be- 
havior of the metric functions practically do not depend 
on the initial value of the exotic scalar field. We have 
investigated also the properties of the mass-function for 
all cases. We want to emphasize the principal difference 
between the collapse of the ordinary scalar field <& with 
e > and the magnetic field H, and the collapse of the 
exotic scalar field \& with eo < and the magnetic field 
H. In the case of the $-field the collapse leads to the for- 
mation of the space-like singularity r = in the T-region, 
but in the case of the \l/-field the singularity r = is a 
time-like and it is in the i?-region. This is correct under 
the assumption which we have made. 

At the end we want to note the following. The consid- 
eration of this paper can be useful also for the analysis 
of the processes inside a BH when it is being irradiated 
by an excotic scalar radiation. 
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fiBr 13 - Br 13 {-A) + Br 13 — A— = 



Formula (|A.3j) gives 



-1 + 4(1 + -) =0, 
From (|A.4j) and (|A.5j) one has 

A-l _ <P_ 
~~A~ ~ ~B 



Formulas (|A.6|) and (|A.7|) give 

A = l, 



5- 



(A.4) 

(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 



Expression (|A.9j) contradicts the conditions d 2 ^ 0. Thus 
the case a — is not possible. 

Now let us consider the case a < 0. To the leading 
order Eqs ©, J3|), q — assume the form 



rr (1) 2 

9tt r 



0, 



,«„(!) 4- „W— =0. 



9tt 9rr + 9r 



(A.10) 



(All) 



APPENDIX A: ANALYSIS OF EQS. ©, © 

We are looking for the solutions to ([3]), (|4]), q = in 
the vicinity r = in the form 



$=Br<>, 
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(i) 



-Ar a , 



(A.I) 



(A.2) 



Let us consider the case a = 0. To the leading order ([3]), 
(U]), q — assume the form 



— A + A + A d B~ r 



2,72 R -l -P-2 _ 



= 



(A.3) 



From (|A.10j) one has 



d 2 

1+ 13=°- 

From (jA.lip and (|A.12[) we obtain 



(A.12) 



(A.13) 



(A.14) 



Expression (j A. 13|) contradicts (|A.14|1 . Thus the case a < 
is not possible. 
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